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ABSTRACT 

In an attempt to constrain evolutionary models of the asymptotic giant branch (AGB) phase at 
the limit of low masses and low metallicities, we have examined the luminosity functions and number 
ratio between AGB and red giant branch (RGB) stars from a sample of resolved galaxies from the 
ACS Nearby Galaxy Survey Treasury (ANGST). This database provides HST optical photometry 
together with maps of completeness, photometric errors, and star formation histories for dozens of 
galaxies within 4 Mpc. We select 12 galaxies characterized by predominantly metal-poor populations 
as indicated by a very steep and blue RGB, and which do not present any indication of recent star 
formation in their color-magnitude diagrams. Thousands of AGB stars brighter than the tip of the 
RGB (TRGB) are present in the sample (between 60 and 400 per galaxy), hence the Poisson noise has 
little impact in our measurements of the AGB/RGB ratio. We model the photometric data with a 
few sets of thermally pulsing AGB (TP- AGB) evolutionary models with different prescriptions for the 
mass loss. This technique allows us to set stringent constraints to the TP- AGB models of low-mass 
metal-poor stars (with M < 1.5 M©, [Fe/H] < —1.0). Indeed, those which satisfactorily reproduce the 
observed AGB/RGB ratios have TP- AGB lifetimes between 1.2 and 1.8 Myr, and finish their nuclear 
burning lives with masses between 0.51 and 0.55 Mq. This is also in good agreement with recent 
observations of white dwarf masses in the M4 old globular cluster. These constraints can be added 
to those already derived from Magellanic Cloud star clusters as important mileposts in the arduous 
process of calibrating AGB evolutionary models. 
Subject headings: stars: general 



1. INTRODUCTION 

The thermally pulsing asymptotic giant branch phase 
(TP-AGB) is both one of the most important and one 
of the more uncertain phases of stellar evolution. Its 
importance resides mainly in its sizeable contribution to 
the integrated light and chemical yields of stellar popula- 
tions, which are essential for the understanding of galaxy 
evolution and the interpretation of the light from dis- 
tant galaxies (e.g . iMaraston et al." '2006'; 'Emi nian et al.l 
[20081 : iConrov et al. 2009; Conroy & Gunn 20To|). Its 
uncertainties derive from a series of circumstances in- 
cluding their complex internal structure, the critical role 
of difficult-to-model processes like convective dredge-up, 
mass loss, circumstellar dust formation and long-period 
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variability, and the scarcity of clear-cut and unequivocal 
observational constraints on their evolution in the most 
immediate universe. 

Gigantic steps are being made in all of these sub- 
jects, but the present situation is that the evolutionary 
timescales of stars in the TP-AGB phase are far from be- 
ing settled. Uncertainties by a factor of a few still exist 
at the extremes of the age-metallicity region allowed for 
AGB stars. This contrasts sharply with other evolution- 
ary phases like the main sequence and red giant branch 
(RGB), for which the evolutionary times are known with 
errors smaller than a few tenths , as demonstrated by a 
large variety of observations (see lGallart et al]|2QQ5l and 
references therein). 

Present day constraints on the TP-AGB evolution are 
largely based on observations of stars in the Magellanic 
Clouds and of the Milky Way (MW). The lifetimes as 
a function of stellar mass, for slightly subsolar metal- 
licities, can be derived for star counts in Magellanic 
Cloud clusters (Frogel et all 119901 : Ivan Loon etahl 120051 : 
iGirardi M arigo 2007al). The same can be done for 
low-mass stars in MW globular clusters. Even for the 
most populous star clusters, however, AGB stars are few 
and their counts are affected by large Poisson fluctua- 
tions; in the case of MW globular clusters, AGB stars 
brighter than the tip of the red giant branch (TRGB), 
and long period variables, practically dis appear from ob- 
served samples at fFe/H] < -1.0 (Frogel EliasI 119881 : 
iFrogel fc Whitelockl[T998[ ). 

Thus, to obtain useful constraints on the AGB lifetime 
from star clusters, it is necessary to either sum the star 
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counts in many clusters in to age and metallicity bins (e.g. 
iGirardi fc MarigQll2QQ7a| ). or embark on a more detailed 
study of the dust, chemical and pulsational properties of 
individual clus ter stars (see e. g. [L ebzelter & Wood 2007; 
iLebzelter et al.n 2008: Kamath et al. 2010 for Magehanic 
Clou d clusters, and L eb zelter et al. 2006; van Loon et al. 
[20061: iMcDonald eFaH ffOOl [20T0I : iBover et al.l l2009bl 
|2010' for MW globulars). Less direct constraints on 
AGB evolution come from integrated cluster proper- 
ties, like their colors and surface brightness fluctuations 
(see e .g. M araston 2005; Pessev et al. 2008; Raimondo 
[2OOI IC^rov et all 120091 : IConrov &: Gunn 2010), which 
by their very nature cannot disentangle different evo- 
lutionary properties like stellar luminosities, lifetimes, 
and chemical types, and hence often provide somewhat 
ambiguous constraints on the numbers of stars at dif- 
ferent evolutionary stages. It may happen, for in- 
stance, that TP-AGB models which do not consider 
the third dredge- up and the formation of ca rbon stars 
(such as the BaSTI ones; ICordier et al.|[2QQ7l ). or mod- 
els in which carbon stars are shifted to the high ef- 
fective te mperature, (Teff) range of 4000 to 3100 K 
(|Conrov et al. 2009; Conrov & Gunn 2010, with their 
logTeff shift of - +0.1 dex for TP-AGB stars at Magel- 
lanic Cloud metallicities), provide acceptable (but spu- 
rious) fits to integrated colors of Magellanic Cloud clus- 
ters. Unfortunately, the properties of these models are 
inconsistent with the most basic observations of over 
10000 C-type AGB stars in the M agellanic Clouds (see 
[Ni kolaev fc Weinberg 20001 iCioni fc Habing 200 3) with 
a significant fraction bein g at 3100 > Teff/K > 2650 
(jGroenewegen et al.ll2009l ). Such high discrepancies and 
gaps in the models remain hidden when only the inte- 
grated properties are examined. 

The problem of low number counts in star clusters 
could, in principle, be circumvented by using AGB star 
counts in entire galaxies, but, in practice, this im- 
plies meeting a series of difficulties. For galaxies with 
moderate-to-high metallicities, the AGB develops at low 
effective temperatures, and consequently near-infrared 
photometry is required to unveil them, as dramatically 
demo n strated by t he M a gellanic Clouds (Frogel et al. 
\Tmi ICioni et al.l [200(1 [Weinberg fc NikolaevI l200lf r 
Moreover, since a no n- negligible fraction of the AGB 
population becomes dust-enshrouded, mid-infrared data 
are necessary f or a compl ete cens us of the AGB 
population (e.g., iBlum et al.P 2006: Bolatto et aP 120071: 
IBover et al.ll2009aD . 

At the limit of very low metallicities, the nearest low 
surface brightness dwarf spheroidal galaxies (dSph) seem 
to be particularly useful for the study of their AGB pop- 
ulations because they appear uncrowded at magnitudes 
accessible with present day 4-m class te lescopes equipped 
with near-i nfrared cameras. Indeed, iGullieuszik et al.l 
(|2008[ ) and iHeld et al.l (|2010D are able to derive impor- 
tant indications from gr ound-based data of the Leo II 
and Leo I dSphs (see also lLagadec et aT]l2008L for Fornax 
and Sagittarius). The main uncertainties in these works 
derive from the small numbers of AGB stars, from the 
high MW foreground contamination, and the uncertain- 
ties in the dSphs s tar formation histories (SFH). Also, 
IBover et al.l (|2009a[ ) call attention to the large fraction 
(30 to 40%) of AGB stars in dwarf irregulars (dirr) which 
are bright in mid-infrared light, but are either missing 



or misclassified in optical studies because they are en- 
shrouded in thick dust shells. It is still not clear whether 
this fraction is also representative of dSphs. 

The possibilities for a better calibration of AGB lu- 
minosities and lifetimes become much wider when milli- 
arcsec resolution im aging is available . A g ood exam- 
ple is provided by iMelbourne et al.l (|2010[ )\ who use 
images taken with the Advanced Camera for Surveys 
(HST/ACS) in the / band (F814W) together with Keck 
Adaptive Optics (AO) in K to study the AGB population 
of the dIrr KKH 98 at a distance of 2.5 Mpc. This work 
clearly demonstrates the utility of moving towards more 
distant resolved dwarf galaxies, for which (1) more stars 
can be observed in a single pointing, even considering the 
small field-of-view available for ACS and AO, and (2) the 
MW foreground becomes dramatically smaller. On the 
other hand, some potential disadvantages of these targets 
are also evident: (a) they present SFHs that are some- 
what more uncertain than the nearest dSphs; (b) there 
is no spectroscopic information to allow a clearcut sep- 
aration between the hottest C- and 0-rich stars (cooler 
stars, instead, are well separated by their infrared col- 
ors), which limits the analysis of aspects related to the 
third dredge- up events in the AGB; (c) the mid-infrared 
photometry is also not available, due to the limited reso- 
lution of Spitzer and AKARI, hence limiting the analysis 
to non-obscured AGB stars; (d) the chance of crowd- 
ing/blending of AGB and RGB stars also increases with 
distance. 

A result common to the low-meta lli city dwarf galaxies 
studie d bv IGullieuszik et al l (|2008f ). IHeld et al.l (|2 01 Of ) 
and Melbourne et al.' ('2010') is that the TP-AGB life- 
tim es from Marigo fc Girardi (200^, hereafter MG07) 
and Marigo et ajj (j2008) models are, apparently, largely 
overestimated. The reason for this is still not clear, 
but is possibly due to the underestimation of the mass- 
loss rates at low metallicities. As a concurring fac- 
tor, the circumstellar extinction presently assumed in 
the models may be somewhat underestimated, as sug- 
gested by the Spitzer observations of nearb y dwarf ir- 
regulars by Jackson et al. (2007a. b) and Bov er et aTl 
(2009a). Regardless, it is likely that the problem is lim- 
ited to low metallicities, since the MG07 TP-AGB mod- 
els have been calibrated in the intermediate-metallicity 
Magellanic Clouds. Their evolutionary behavior at very 
low metallicities can be considered as either educated 
guesses, or extrapolations of the behavior met at interme- 
diate metallicities, since they are a result of the straight 
application of theoretically-uncertain metallicity depen- 
dences. 

In this paper, we use HST optical data from the 
ANGST survey to derive constraints on the AGB evolu- 
tion for a set of metal-poor galaxies in which the recent 
star formation activity is generally very low, if not ab- 
sent. Our target is to derive constraints on the optically- 
visible TP-AGB lifetimes of low-mass, metal-poor stars. 
The advantage of using a set of Mpc-distant galaxies is 
clear: we will have thousands of AGB stars in our sam- 
ples, reducing Poisso n noise to a minimum. The price to 
pay, as in the case of lMelbourne et al.l ()2010[ ). is that our 
analysis can be affected by uncertainties in the SFHs of 
individual galaxies. The hope however is that these SFH 
errors become less relevant when averaged over a large 
sample of galaxies. 
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In Sect. [2] we describe the data used in this work. In 
Sect. [3] we model the data using two different sets of 
TP- AGB models, deriving clear quantitative constraints 
on them. We then provide a revised set of TP- AGB 
tracks that brings models and observations into better 
agreement. 

2. THE DATA 

2.1. The ANGST/ANGRRR photometry database 

We used data from the ACS Nearby Gala xy Survey 
Treasury (ANGSTJB (jPalcanton et al.il2QQ9[ ) and from 
the "Archive of Nearby Galaxies: Reduce, Reuse, Recy- 
cle" (ANGRRRjE databases. These two programs have 
archived stellar photometry for tens of millions of stars 
in nearby galaxies outside the Local Group, based on im- 
ages taken with the ACS and the Wide Field Planetary 
Camera 2 (WFPC2) on the HST. The imaging primarily 
used blue observations in the F475W, F555W, or F606W 
filters, and F814W for the red filter. The photome- 
try was performed using DOLPHOT0 and HSTPHOT0 
(for ACS and WFPC2 im aging, respective l y; iDolphinI 
2QQQD. as descr i bed fu lly in iDalcanton et al.l (j2QQ9[ ) and 
Williams et al.l (|2QQ9f ) . The depth of the resulting color- 
magnitude diagrams (CMD) varies with distance, crowd- 
ing, and exposure time, but reaches several magnitudes 
below the TRGB in all cases, and below the red clump 
(M/ —0.5 mag) in all but 2 of the cases we analyze 
here. 

For this paper, the photometry from the first data re- 
lease for ANGST and ANGRRR has been supplemented 
with extensive artificial star tests. For each of the fields 
studied, a minimum of 100,000 artificial star experiments 
were performed to determine the completeness and er- 
rors as a function of color and magnitude. In each test, 
one star of known color and magnitude was added to 
the data, and the photometry routine was re-executed to 
determine the difference between the input and output 
magnitudes if the star was recovered. 

SFHs were determined by fitting the observed CMDs 
with the MATCH fitting package (jPolphinl 120021 ) . The 
parameters of MATC H were set to match those us ed by 
iWilhams et HI (^2009"). In brief, a "Salpeter (1955) IMF 
was used to generate a base set of Hess diagrams us- 
ing the i sochrones of IGirardi et al.l (j2002, with updates 
in M arigo et al.l l2008f ) and the error and completeness 
from our artificial star tests. The best-fit combination 
of this base set of diagrams was then dete rmined us- 
ing the statistics described in Dol phinI ()2002l ). First, the 
fit was performed on the full data set brightward of the 
50 % completeness limit as determined by our artificial 
star tests. Then, the fit was repeated excluding all data 
brightward of the TRGB. This second fit allowed us to 
quantify the effects of the bright AGB stars on the re- 
sulting SFH (as discussed in Sect. 13. 2p . 

We further analyzed the CMDs to estimate extinctions 
Ay^ distance moduli (m — M)o, and the magnitude of 
the TRGB. The total extinction and distance moduli are 
determined automatically by MATCH, based on the op- 

http:/ / www.nearbvgalaxies.Qrg, 

'http: / / archive.stsci.edu / prepds / angst / 

http:/ /archive. stsci.edu/prepds/angrrr 
^'^ http://purcell.as.arizona.edu/dolphot 
http:/ /purcell. as. arizona.edu/hstphot 
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Fig. 1. — Top panel: F814W-band bolometric corrections as 
a function of Teff, for both [Fe/H] = (continuous line) and 
[Fe/H] = —2 (dashed line). Bottom panel: The Teff-color re- 
lations for both F475W-F814W and F606- F814W, for the sam e 
metallicities. These relations are described in lGirardi et al.l (|2008I V 
The small kinks at Te g ^ 3800 K correspond to th e transition 
between Castelli &; Kumci (l2003l . for log g = 2) and IFluks et al.1 
(1994) model atmospheres. 

timal reddening and distance needed to reproduce the 
obser ved CMD (see the full discussion in Williams e t al] 
assuming Ry = 3.1 at optical wavelengths. Un- 
certainties are determined by identifying the range over 
which Av and (m — M)o can vary without producing a 
statistically significant reduction in the quality of the fit 
to the observed CMD. 

2.2. The galaxy sample and their SFHs 

We chose galaxies for our sample based on two primary 
criteria - metallicity and mean stellar age. The first re- 
striction results from our desire to analyze only the galax- 
ies in which the ANGST/ANGRRR optical data contains 
a nearly complete census of the AGB population. This 
requirement naturally limits our work to galaxies with 
metal-poor populations (say [Fe/H] < — 1.2). At these 
metalHcities, the upper AGB develops at Teff > 4000 K 
and colors F475W-F814W< 2.4, or F606W-F814W< 
1.2, and still falls at the plateau of the BC-Teff relations 
(Fig.[T]). Bolometric corrections for low- metallicity AGB 
stars are thus close to their minimum values, ensuring 
that the brighter stars really do appear at the smaller 
magnitudes. In contrast, for redder, high- metallicity 
AGB stars, bolometric corrections become significant 
and increase steadily with color. Redder AGB stars thus 
become progressively less accessible in the optical, and 
instead require near-infrared photometry to extend the 
complete AGB counts towards the region of lower Teff, 
which corresponds to more metal-rich galaxies. 

We further restrict the galaxy sample to regions that 
have a simple-to-interpret SFH, with AGB stars limited 
to a well-defined region of the age-metallicity plane. This 
request limits us to galaxies dominated by old stars, for 
which there is no sign of recent or intermediate-age star 
formation. Thanks to the steep mass-main sequence life- 
time relation, the absence of populations younger than 
3 Gyr ensures that the evolved stars have initial masses 
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Fig. 2. — The data and models used in this work, for all galaxies and galaxy regions in our sample. Left panels: The observed CMD 
from ANGST/ ANGRRR. The horizontal lines mark the TRGB and upper and lower magnitude limits of the stars considered to be in the 
upper RGB and AGE. Their total numbers are marked inside the boxes. Middle panels: The same as in the left panel, but for the model 
CMD from TRILEGAL, with the default set of TP-AGB models (Marigo et al. 2008) and scaled to present the same number of upper 
RGB stars. Right panels: A comparison between the observed (black) and model (red or grey) luminosity functions. 



confined to a narrow interval, of rougiily 0.8 to 1.4 Mq. 
The AGB star counts in tiiese galaxies should then pro- 
vide clear constraints on the evolution of low-mass stars. 

Based on these two conditions, we have selected a sam- 
ple of galaxies for which a visual inspection of the CMD 
revealed: (1) an almost vertical RGB, with the TRGB at 
F475W-F814W< 2, or F606W-F814W< 1 so as to in- 
dicate low metallicity and BCf814W ^ 0.5 (Fig. [T]), and 
low foreground extinction, and (2) no evidence for re- 
cent and intermediate-age star formation, as indicated 



by an absence of any younger main sequence and helium 
burning sequence above the red clump/HB. This latter 
restriction eliminates galaxies with star formation in the 
most recent 0.5 Gyr, but may admit galaxies with some 
amount of star formation at intermediate ages. A full 
analysis of the SFHs can be found in Weisz et al. (in 
preparation). 

Although our selection considered both ACS and 
WFPC2 data, it turns out that the final sample contains 
ACS data only. 
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Fig. E] continued. — 

The left panels in Fig. [2] illustrate the CMDs for ah the 
selected galaxies and galaxy regions. Their basic prop- 
erties and parameters are listed in Table [TJ 

2.3. Selecting AGB and RGB stars 

To identify the location of RGB stars , we adopted the 
magnitude of the TRGB (mTRGB) from lDalcanton et al.l 
|2QQ9i) . who identified the TRGB in uncrowded low ex- 
tinction regi ons of each g alaxy using the edge-detection 
method of Mendez et al. (2002). The TRGB magni- 
tude can also be converted into a distance mod ulus, us- 
ing th e foreground extinction derived by Schl egel et al.l 
(|1998l ). and t he absolute ni a gnitu de of the TRGB in the 
isochrones of iMarigo et al.l (|2008f ) at the observed color 
of the RGB stars used to derive the TRGB. The TRGB- 
based distance modulus does not necessarily agree with 
the one reported by MATCH, since the latter can be bi- 
ased towards larger distances in an attempt to better fit 
the magnitude of the well -populated red clum p, which is 
known to be faint in the Marigo et al.l (|2QQ8f ) isochrones 
currently used in MATCH. In what follows, we use the 
(m — M)o and Ay derived from MATCH for generating 
artificial CMDs. However, for isolating AGB stars, we 
use the empirically determined value of mxRGB , as given 
in Table E 

An interval of two magnitudes above the TRGB defines 
the "AGB sample". In all cases, it includes the bulk of 
(if not all) stars brighter than the TRGB. Their number, 
^AGB, is typically between 60 and 400 per galaxy. 

The RGB sample, instead, is defined over two mag- 
nitudes below the TRGB. This sample has, typically, a 
completeness above 95 %, and about 30 times more stars 
than in the AGB sample. 

The number ratios between AGB and RGB stars, 
^agb/^rgb, are presented in Table [H together with 
the Icr random errors. 

The AGB sample as above defined is contaminated by 
a few RGB stars at its faintest magnitude bins, because 
of the scattering of stars to brighter magnitudes by pho- 
tometric errors and binaries. For the galaxies with the 
smallest A^agb/^rgb, the effect is such that this ratio 
can be increased by up to 20 % with respect to its true 
value. As we will see below, this is a marginal effect con- 
sidering the large discrepancies - of a few times - between 
present models and the observations. Moreover, the scat- 



tering of RGB stars to brighter magnitudes is fully taken 
into account in our simulations (Sect. [3]), so that no in- 
consistency results from counting a few RGB stars in 
the AGB sample. On the other hand, the RGB sample 
contains both the initial sub-luminous section of the TP- 
AGB phase and the excursions to low-luminosities driven 
by thermal pulses, as well as stars leaving the early- AGB 
phase. Based on evolutionary models, we estimate that 
this contamination, in the worst cases, is just a few per 
cent. Moreover, the numbers of early- AGB stars are ex- 
pected to be quite insensitive from the uncertainties in 
the mass-loss prescriptions, which instead plague the TP- 
AGB stars above the TRGB. 

2.4. Uncertainties in the SFH and AGB/RGB ratio 

The measured A^agb/^rgb ratios have values com- 
prised between 0.023 and 0.050. Random errors are, in 
all cases, smaller than 20 %. The mean value derived by 
adding all stars in the sample is 0.0344 ±0.0007. Clearly, 
the random noise in the relative number of AGB stars 
is less of a concern in our dataset. The main point of 
concern, instead is in the correctness of the SFHs for the 
host galaxies, which is crucial for the interpretation of 
the A^agb/A^rgb ratios. 

As we show below, the current isochrones overpopulate 
the AGB sequences, potentially biasing our SFHs to low 
star formation rates and intermediate ages when these 
stars are included in the fit of the CMD. We therefore re- 
derive the SFHs excluding AGB stars, taking care to keep 
all other aspects of the fitting identical. Specifically, we 
exclude CMD regions above the TRGB, and those with 
either F606W-F814W> 3.5 or F475W-F814W> 4.5 
(depending on the filters used in the observations) from 
our fitting. The resulting SFHs correspond to the "no- 
AGB" cases in Table El They differ little from the default 
"with-AGB" cases, indicating that the structure of the 
red clump has far more influence on the inferred SFH at 
intermediate ages with the weighting scheme currently 
employed by MATCH. 

We also have applied the "Z-inc" option to derive SFH. 
It consists of limiting the number of free parameters in 
our fits by forcing the fit to only attempt solutions where 
the metallicity remains constant or increases with time. 
The history returning the best fit to the CMD is then 
selected. This process is important in some cases where 
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TABLE 1 

Basic parameters of our galaxy sample. 



Target / name 


Filters 


Ay 




''''50%complete 


(m — M)n 


R B 




Nag-b 

-^RGB 


ESO294-010 


F606W,F814W 


0.018 


22.397±0.008 


28.25 


26.434 


1950 


66 


0.034±0.004 


UUKJ LLo 


r 4 / W ,r oi4W 


U.Ubo 


OQ QQ^y -Un r\oci 




z / .o4y 


z (UZ 




U.UoOzLU.UU4 


DD044 


F475W,F814W 


0.129 


23.499±0.015 


27.61 


27.454 


4966 


205 


0.041±0.003 


ESO540-032 


F606W,F814W 


0.064 


23.733±0.020 


27.83 


27.743 


3627 


117 


0.032±0.003 


M81F6D1 


F606W,F814W 


0.241 


23.855±0.010 


27.85 


27.737 


2326 


60 


0.026±0.003 


M81K61 


F606W,F814W 


0.226 


23.823±0.042 


28.05 


27.716 


8181 


263 


0.032±0.002 


IKN 


F606W,F814W 


0.181 


23.869±0.019 


26.97 


27.786 


7854 


184 


0.023±0.002 


DD078 


F475W,F814W 


0.066 


23.837±0.017 


27.55 


27.818 


9559 


394 


0.041±0.002 


SCL-DEl 


F606W,F814W 


0.046 


24.116±0.023 


28.41 


28.110 


1922 


96 


0.050±0.005 


DD071 


F606W,F814W 


0.303 


23.881±0.019 


28.05 


27.740 


6540 


255 


0.039±0.002 


M81F12D1 


F606W,F814W 


0.442 


23.982±0.013 


28.03 


27.749 


7294 


236 


0.032±0.002 


M81K64 


F606W,F814W 


0.165 


23.907±0.009 


28.42 


27.852 


6876 


291 


0.042±0.003 


total 












65771 


2265 


0.0344±0.0007 



Note. — Ay and (m — M)o come from the best MATCH solution. 



the photometry depth or the number of stars are insuffi- 
cient to reliably separate the effects of age from those of 
metallicity. 

The initial columns in Table [2] summarize the results 
in terms of SFH, presenting the fraction of the SFH in 
the age intervals from to 1 Gyr, and from 1 to 3 Gyr, 
for all cases of SFH recovery that were tested. It turns 
out, as expected, that all galaxies but SCL-DEl present 
a very small amount of their SFH at ages < 1 Gyr, and 
roughly 10 % at 1 - 3 Gyr. The SFHs for the different 
cases (with and without AGB stars, default or Z-inc) 
do agree well considering the typical error bars in this 
kind of measurement, again with the exception of SCL- 
DEl for which the Z-inc cases provide significantly older 
SFHs. 

3. MODELING THE DATA: METHOD AND RESULTS 
3.1. Code and method 

To model the ANGST/ANGRRR data, we use a recent 
version of the TRILEGAL code (IGirardi et al.l [2005). 
which generates multi-band mock catalogs of resolved 
stellar populations for a given distribution of distances 
and extinctions, following some specified star formation 
history and age-metallicity relation. The original code 
has been expanded in several aspect s to deal with the 
compl ex features of TP-AGB stars (IGirardi & Marigd 
I2007U and later work), including: different bolometric 
corrections and Teff-color relations for 0-rich and C-rich 
stars, luminosity (L) and Teff variations driven by ther- 
mal pulses, and obscuration by circumstellar dust. In 
order to take these effects onto account, the code keeps 
track of a series of stellar parameters, like the surface 
chemical composition, mass loss, and the period and 
mode of the long-period variability. The dust composi- 
tion has been assumed to be the 60 % silicate plus 40 % 
AlOx for 0-rich stars, and 85 % amorphous c arbon plus 
15 % SiC for C-rich stars from iGroenewege 3(^006). 

The ste llar evolutiona r y trac ks are the same ones con- 
tained in iMarigo et al.l (|2008l) isochrones and used by 
MATCH to derive the SFH of our galaxy sample. The 
transformations from L and Teff t o the H ST photome- 
try are described in IGirardi et al.l (|2008f ): however, we 
have updated it t o use the late s t tran sformations for 
C-type stars from lAringer et al.l ()2009f ). and the total 
ACS/WFC throughput curves and zeropoints appropri- 
ate for post- July 2006 observations (Mack et al.. .2007: 



iBohiidlmrliFI . 

Alternative sets of TP-AGB tracks have been specifi- 
cally calculated for the present work, as described below 
in Sect. 13.41 TRILEGAL allows these tracks to be re- 
placed quickly and with minimal human effort. 

We first start modelling the data for the complete sam- 
ple of galaxies in Table [H using the same stellar models 
which were already used to derive the SFH, but using 
TRILEGAL instead of MATCH. TRILEGAL simulates 
the photometry starting from the SFH file, and shifts the 
data to the right distance modulus and extinction. The 
simulations include stars up to 2 mags fainter than the 
faintest observed one. 

In order to properly account for the real completeness 
and photometric errors in our simulations, we proceed 
as follows. For each simulated star, an artificial star 
of similar color and magnitude is randomly extracted 
from the existing catalog of artificial stars from HST- 
PHOT/DOLPHOT. If that artificial star has been de- 
tected by the photometry software, the differences be- 
tween the input colors and magnitudes, and the output 
ones, are applied to the simulated star, otherwise the 
same object is thrown away. 

3.2. Results using MG07 TP-AGB tracks 

Results from this exercise are shown in Fig. [2l Notice 
that the simulations are forced to have a number of bright 
RGB stars (within 2 mag of the TRGB) consistent with 
the observed one within la. One can notice that these 
simulations resemble very much the observations, except 
for the TP-AGB regime, for which there is a clear excess 
of simulated stars, by factors that can be as large as 6 or 
7. The numerical results are tabulated in Table [2] below. 
It is also evident that the simulated AGB stars reach 
much brighter magnitudes than the observed ones. The 
excess of AGB stars is^probably relate d to w h at has been 
already detecte d bv iGu llieuszik et alj (|2008[ ) , iHeld et all 
(2010), and M elbourne et al. (2010). 

We have checked that these results depend little on the 
use that MATCH makes of AGB stars. Indeed, as shown 
in Table El the SFH derived when one hides the AGB 
stars in MATCH are very similar to those in which they 
are included, and do not produce great changes in the 
predicted AGB/RGB ratios. Moreover, the comparison 

■"^^ http:/ /www. stsci.edu/hst/acs/analysis/zeropoints 
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TABLE 2 
The SFR and AGB/RGB ratio 



Target / name 



Which SFH 



A^agb/A^rgb 





AGB? Zinc? 


-'^< IGyr 
SFRtot 


-t^l — 3Gyr 
SFRtot 


{[Fe/H]> 


observed 


Ma08 


case A 


case B 




Y 


N 


0.01 


0.07 


-1.68 







149 





020 





020 




Y 


Y 


0.01 


0.07 


-1.63 







131 





021 





019 




N 


N 


0.02 


0.08 


-1.75 







133 





025 





019 




N 


Y 


0.02 


0.09 


-1.74 







155 





021 





022 


DD01 1 S 


Y 


N 


0.02 


0.23 


-1.55 







147 





068 





073 




Y 


Y 


0.01 


0.17 


-1.57 







144 





070 





062 




N 


N 


0.03 


0.16 


-1.44 







162 





076 





068 




N 


Y 


0.03 


0.17 


-1.42 







155 





066 





075 


DD044 


Y 


N 


0.01 


0.21 


-1.20 


n 041 _Ln nno 





147 





070 





069 




Y 


Y 


0.01 


0.19 


-1.13 







131 





039 





050 




N 


N 


0.02 


0.11 


-1.13 







150 





064 





067 




N 


Y 


0.02 


0.21 


-1.14 







146 





065 





066 




Y 


N 


0.03 


0.05 


-1.42 


n nooin nno 
yj . L/«jZ/_nj . yjyj^j 





146 





055 





049 




Y 


Y 


0.02 


0.03 


-1.47 







166 





058 





045 




N 


N 


0.03 


0.07 


-1.56 







162 





043 





042 




N 


Y 


0.02 


0.03 


-1.48 







162 





044 





046 




Y 


N 


0.01 


0.04 


-1.12 


n 026+0 00s 

yj . I f Jill ,\JyjKJ 





133 





024 





026 




Y 


Y 


0.01 


0.02 


-1.25 







145 





033 





030 




N 


N 


0.02 


0.06 


-1.17 







137 





030 





028 




N 


Y 


0.02 


0.01 


-1.20 







149 





025 





026 


MSI Kfil 


Y 


N 


0.01 


0.04 


-1.12 


n nooin nno 





112 





020 





021 




Y 


Y 


0.01 


0.02 


-1.25 







122 





034 





035 




N 


N 


0.02 


0.06 


-1.17 







135 





039 





043 




N 


Y 


0.02 


0.01 


-1.20 







147 





044 





047 


IKN 


Y 


N 


0.02 


0.03 


-1.11 


n nooin nno 
yj . yjZiKj -\—yj . yjyjZi 





116 





029 





031 




Y 


Y 


0.02 


0.05 


-1.49 







124 





028 





031 




N 


N 


0.03 


0.00 


-1.17 







160 





046 





040 




N 


Y 


0.04 


0.00 


-1.55 







158 





033 





033 


DD078 


Y 


N 


0.01 


0.14 


-1.12 


n n4-| _Ln nno 





127 





069 





071 




Y 


Y 


0.01 


0.11 


-1.10 







156 





089 





088 




N 


N 


0.01 


0.16 


-1.17 







150 





081 





079 




N 


Y 


0.01 


0.15 


-1.09 







164 





095 





102 




Y 


N 


0.10 


0.27 


-0.79 


n nirnin nnc; 
yj . yj'jyj -\—yj . yjyjo 





314 





269 





251 




Y 


Y 


0.06 


0.11 


-1.57 







358 





274 





275 




N 


N 


0.06 


0.38 


-0.98 







136 





061 





052 




N 


Y 


0.08 


0.05 


-1.63 







212 





105 





099 


DD071 


Y 


N 


0.01 


0.12 


-1.05 


n noQ_|_n nno 





121 





021 





023 




Y 


Y 


0.02 


0.08 


-1.19 







115 





029 





030 




N 


N 


0.02 


0.11 


-1.06 







141 





040 





046 




N 


Y 


0.01 


0.07 


-1.17 







135 





038 





043 


M81F12D1 


Y 


N 


0.01 


0.08 


-1.17 


0.032±0.002 





127 





020 





018 




Y 


Y 


0.01 


0.03 


-1.29 







125 





023 





024 




N 


N 


0.01 


0.05 


-1.15 







142 





034 





032 




N 


Y 


0.01 


0.02 


-1.19 







130 





023 





024 


M81K64 


Y 


N 


0.01 


0.08 


-1.25 


0.042±0.003 





134 





022 





023 




Y 


Y 


0.01 


0.08 


-1.32 







138 





030 





025 




N 


N 


0.01 


0.07 


-1.26 







143 





038 





033 




N 


Y 


0.01 


0.06 


-1.30 







146 





035 





039 








Results for mid-to-lowest density rings 


of galaxies 














ESO540-032 


N 


N 


0.03 


0.06 


-1.49 


0.028±0.004 





162 





059 





058 


DD0113 


N 


N 


0.02 


0.13 


-1.45 


0.037±0.005 





153 





071 





052 


M81K61 


N 


N 


0.01 


0.11 


-1.21 


0.029±0.003 





143 





043 





042 


DD078 


N 


N 


0.02 


0.12 


-1.24 


0.038±0.003 





169 





070 





076 



with models in which the metallicity is forced to not de- 
crease with the galaxy age (the Z-inc option in MATCH) 
also produces similar results. We can only conclude that 
the problem resides in the AGB stellar models, and not 
in the observed samples or process of SFH recovery. 

The observed galaxies always present, overall, a small 
fraction of intermediate age stars (Table [2]), or some 
moderately-high metallicities, so that one might think 
that the results do not actually correspond to old metal- 
poor populations as it was intended to be. We know 
however that in dwarf galaxies the youngest (and more 
metal-rich) star formatio n is, as a rule, concentrated in 
the g alaxy centers (e.g. Weisz et"ani2QQ8l : iStinson et al.l 
l2QQ9l and references therein; also Gilbert et al. in prepa- 



ration). Therefore, for a subsample of our dwarf galaxies, 
we have selected the external rings for which the SFH is 
expected to be older than for the overall galaxy. The 
results are amended at the end of Table [21 for the no- 
AGB no-Z-inc option only. Although these outer regions 
are indeed slightly older than the entire galaxies, they 
show essentially identical A^agb/^rgb ratios as the in- 
ner regions, suggesting that old populations dominate 
the CMDs at all radii. Therefore, the basic result for the 
excess of predicted AGB stars, by a factor of about 4, 
still remains in these data. 



3.3. What is the problem with current TP- AGB models? 
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The above-mentioned results were obtained with the 
use of MG07 TP-AGB evolutionary tracks, which consti- 
tute a quite non-standard grid of such models. They have 
considered, for the first time, features like the changes in 
molecular opacities and mass-loss rates as the surface 
composition of TP-AGB stars passes from O- to C-rich, 
and the increase in the mass-loss rates as the long-period 
pulsation switches from the first overtone to the funda- 
mental mode. Added to these improved prescriptions, 
there were attempts to calibrate the poorly known pa- 
rameters of the models via the fitting of observational 
data. To be explicitly considered in the fitting were (1) 
the lifetimes of AGE stars as a function of stellar mass 
as derived from star counts in Magellanic Cloud clusters, 
and (2) the C-star luminosity functions in both Clouds. 
After this calibration, other quantities were checked, as 
for instance the integrated colors of Magellanic Cloud 
clusters, and the initial-final mass relation derived from 
white dwarfs in the solar neighborhood. All these com- 
parisons revealed an overall improvement in fits to data, 
even if there remained some problems, like for instance 
the integrated V — K and J — K colors being apparently 
too red as compared to the cluster data in the a ge interva l 
from 10^ to 4 X 10^ yr (see the appendix in iMarigo et al.l 

[lool). 

One can separate the prescriptions used to build MG07 
models into two broad groups: those which depend on 
the composition of stellar atmospheres - and especially 
on the C/0 ratio - and those which do not. The first set 
of prescriptions is probably of secondary relevance for the 
present work, since we are dealing with stars in metal- 
poor and old stellar systems, most of which do not suffer 
convective dredge-up on the TP-AGB and hence are ex- 
pected to be predominantly 0-rich (of spectral types K 
and M). For the stars we are considering here, the TP- 
AGB lifetime and termination luminosity are crucially 
determined by the process of mass loss. 

For 0-rich stars, MG07 adopte d a luass-lo ss prescrip- 
tion heavily based on' Bowen fc WillsonI (1991) dynamical 
model atmospheres for fundamental- mode pulsators, and 
with a dependence on metallicity as derived from Willson 
(2000). The initial period of TP-AGB mass-loss as first- 
overtone pulsators was based on a set o f relati ons de- 
rived from very few exploratory models bv' BowenI p988.) . 
The procedure however was not straightforward, involv- 
ing a large number of assumptions and fitting relations, 
caused by the different assumptions and ranges of param- 
eters (e.g., the different underlying radius-luminosity- 
mass relation) between the original models from Bowen 
and Willson, and the synthetic TP-AGB models. A de- 
tailed review of the literature, together with the results 
of Fig. [21 suggests a complete revision of this scheme, as 
we now describe. 

3.4. New TP-AGB tracks 

Our main problem is to reduce the lifetimes of TP- 
AGB tracks which, due to their low mass and metallic- 
ity, present relatively low luminosities and hot temper- 
atures. Under such conditions, the classical dust-driven 
winds typically found in more luminous and cool AGB 
stars might be less efficient. In any case, we need to 
introduce some description for the "pre-dusty" winds 
(with rates Mpre-dust), i-e. before radiation pressure on 



dust grains becomes th e main driving agent of mass loss 
(jElitzur fc Ivezicll2001f ). Among the few such formulas 
which have been proposed in t he literature, we dec i ded to 
test the semi-empirical one bv lSchroder fc Cuntzl (l2005h 



LR 



pre— dust — V 



teff 



M V4000K 



3.5 



9Q 

4300 



[Moyr- 



(1) 

where i?, M, and L are the stellar radius, mass, and 
luminosity expressed in solar units; g and Qq indicate 
the stellar and sol ar surface grav ity, respectively. This 
formula is a kind of lReimersI (|1975[ ) law and was originally 
derived to describe the mass loss suffered by red giants 
under the assumption that the stellar wind originates 
from magneto-acoustic waves operating below the stellar 
chromosphere. The fitting pa rameter 77 is set to 0.8 x 
10-^^, a value calibrated by ISchroder fc Cuntzl (j2005D 
to reproduce the morphology of horizontal branches in 

globul ar clusters. 

The ISchroder fc Cuntzl (|2005 ^ relation is employed 
from the beginning of the TP-AGB phase until a crit- 
ical stage is met, i.e., the attainment of the minimum 
mass- loss rate, M^^^, required for the development of a 
dust-driven wind. At each time step M"^^^^ is evaluated 
numerically following the analysis by iGail fc Sedlmavrl 
(1987) , from the condition 



L K 



a 



AttcGM 



1. 



(2) 



which expresses the balance between the outward force 
caused by radiation pressure on dust and the inward 
gravitational pull of the star. Here c is the light speed, G 
is the gravitational constant, and is the flux- averaged 
mass extinction coefficient of the gas-dust mixture. Fol- 
lowing Ferrarotti fc Gaill (|2006t) it is reasonable to as- 
sume that in the dust-driven outflow 



^gas 



J2 



ust,i 



(3) 



where the opacity contributions from the gas, A^gas, and 
dust species, K,dust,i^ are expressed as Rosseland means, 
while fi represents the condensation degree of the i^^ 
dust species under consideration. 

In our calculations we consider different dust types 
depending on the star chemical type: silicates (pyrox- 
ene, olivine and quartz) and iron in the case of M-stars 
(C/0 < 1), iron in the case of S-stars (C/0 ~ 1), 
and carbon, silicon carbide, and iron in the case of 
C-stars (C/0 > 1). The dust op acities are evalu- 
ated th rough the relations provi ded bv IGail fc Sedlmavr 
(1999), iFerrarotti fc Gaill (|200lD . and iFerrarotti fc Gail 
(p002). and the corresponding condensation degrees are 
estima ted with the analytic fits proposed by Ferra rottil 
([2003h . 

These latter depend on the current M, so that the 
minimum mass- loss rate of a dust-driven wind, M^^^, is 
found iteratively through Eq. (|2]). 

As soon as ikfpre— dust ^ ^dust' star enters the dust- 
driven wind regime, which is treated according to two 
formalisms: 



case A: based on iBowen fc Wills on (1991), but 
relaxing the metallicity dependence suggested by 
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iWillsonl (|2QQ(I : 

• case B: based on "Bedijn' (1988'), but with a some- 
what different calibration of the parameters. 

Case A for dust-driven winds is the same as in MG07, 
but for the correction factor describing the explicit metal- 
licity dependence, which is left out in the new TP- AGB 
models. Nonetheless, an intrinsic metallicity effect still 
remains in the computed mass-loss rates, via the stellar 
surface parameters R and Teff. In fact at lower metallici- 
ties the atmospheres of AGB stars tend to be hotter and 
more compact. In practice, case A is meant to explore 
the hypothesis that the dust-driven winds may depend 
only mildly on the initial metal content. 

Case B clo sely resembles the approach developed by 
iBedijnl (|1988f ). to which the reader is referred for all de- 
tails. Briefly, assuming that the wind mechanism is the 
combined effect of two processes, i.e., radial pulsation 
and radiation pressure o n the d ust grains in the outer- 
most atmospheric lavers. lBedijnl (jl988) derived a formal- 
ism for the mass- loss rate as a function of basic stellar pa- 
rameters, M, R, Teff, and the photospheric density pph- 
Similarly to iBedijnl (j 19881 see his Fig. 1 and appendix 
B) the free parameters have been calibrated on a sample 
of Galactic long-period variables with known mass-loss 
rate, pulsation period, stellar mass, radius, and effective 
temperature. More details about the fit procedure will 
be given elsewhere (Marigo et al., in preparation). 

It should be noted that, in our calculations, i?, Teff, and 
Pph are derived from numerical integrations of complete 
envelope models extending from the photosphere down 
to the degenerate C-0 core (see Marigo et al. 1999, for 
details). A key prerogative of our TP- AGB code is that 
at each time step, low- temperature opacities are com- 
pute d, for the first time, on- t he-fly with the iESOPUS 
tool ([Marigo fc Aringerll2QQ9f ). thus assuring a full con- 
sistency with the surface chemical composition. In this 
way we avoid the loss in accuracy that otherwise must 
be paid when interpolating on pre-computed opacity ta- 
bles (the standard approach in stellar evolution mod- 
els). This is important since low-temperature opacities 
are crucial in determining the position of a giant in the 
Hertszprung-Russell diagram. The opacities of carbon- 
rich stars (C/0 > 1), for instance, may be significantly 
higher, on average, than in oxygen-rich stars (C/0 < 1), 
so that the atmospheres of carbon-rich stars are usually 
less dens e (and mor e extended) than those of oxygen-rich 
stars (.Marigoll2QQ2h . 

Figure [3] compares the MG07 results with the new ones, 
in terms of (i) duration of the TP- AGB phase and (ii) 
initial-final mass relation. It is evident the significant 
reduction (by a typical factor of 3 — 4) of the TP- AGB 
lifetimes for initial masses Minitiai ^ 0.9 M©, as a conse- 
quence of the new prescriptions for the mass- loss rates. 
The shortening of the TP-AGB lifetimes derives from 
two concurring factors, i.e.: the non- negligible mass- loss 
rates already prior the onset of the dust-driven wind, and 
the higher efficiency during the dusty regime compared 
to MG07. 

This fact is illustrated in Fig HI showing the predicted 
evolution of the mass-loss rate for a model with ini- 
tial mass Minit = 0.832 M© and metallicity Z = 0.001, 
according to the three different choices of mass-loss 




0.6 0.8 1 1.2 1.4 1.6 

Fig. 3. — Comparison of AGB models with initial metallicity Z = 
0.001 and for three different mass-loss prescriptions, as described 
in the text. Top panel: TP-AGB lifetimes as a function of the 
initial stellar mass. Bottom panel: Final core masses left after the 
ejection of the envelope as a function of the initial stellar mass. 
The hatched rectangle shows the initial-final mass relation for the 
Galactic globular cust er M4, a ccording to the recent white dwarf 
measurements by Kalir ai et al.l f2009). 

formalisms discussed here. The evolutionary time is 
counted since the first thermal pulse up to the ejection of 
almost the entire envelope (calculations are stopped as 
soon as the envelope mass falls below 0.01 Mq). One crit- 
ical feature common to cases A and B is that, despite the 
low metallicity, the mass- loss rate Mpre-dust predicted by 
the lSchroder fc Cuntd (|2005l ) relation is already efficient, 
quickly increasing from ^ 10~^ to 5 x lO~^M0yr~^ 
in a few thermal pulses. These rates are high enough to 
determine the ejection of the entire stellar mantle even 
before the development of the dust-driven wind (case B; 
bottom panel), or to favor its earlier onset (case A; mid- 
dle panel). In general the use of mass- loss prescription A 
leads to somewhat shorter TP-AGB lifetimes than case 
B. 

Consequently, compared to MG07, the models with 
the new mass loss descriptions predict a sizable reduc- 
tion of the final masses, that are now mostly com- 
prised in the interval 0.51 < (Mfinai/M©) < 0.55 for 
0.70 < (Minitiai/Mo) < 0.90 (bottom panel of Fig. [3]). 
The most important result is that with no attempt to 
tune the mass loss, the new TP-AGB models are able 
to naturally recover the empirical initial-final mass re- 
lation of p opulation II stars (b ottom panel of Fig. [3]) as 
derived bv lKalirai et al.l (|2009[ ) from the first direct mass 
measurements of individual white dwarfs in the Galac- 
tic globular cluster M4. We have also verified that all 
metallicity sets in the range 0.0001 < Z < 0.001 intersect 
the empirical data. 

At larger initial masses, 0.9 < Minitiai /^o ^ 1-5, in- 
stead, there is a general agreement between the three 
sets of models, which are all predicted to make the tran- 
sition to the C-rich domain (C/0 > 1), given the higher 
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efficiency of the third dredge-up at the low metallicities 
here considered. In this case the short TP-AGB lifetimes 
are controlled by the higher mass-loss rates for C-stars 
(following the reduction of Teff as soon as C/0 exceeds 
unity), while the flattening of the initial- final mass re- 
lation is mainly shaped by the deep dredge- up events. 
Interestingly, the three mass-loss prescriptions adopted 
here converge to similar results. 

A general finding is that low-mass low-metallicity mod- 
els (Minitiai ^ I.OMq) would remain 0-rich for all their 
TP-AGB evolution, experiencing non-negligible mass 
loss already before reaching the conditions to activate a 
dust-driven wind (not excluding that some passive dust 
could actually be present), so that the mass- loss mech- 
anism on the TP-AGB should be ascribed to a different 
(unknown) driver, e.g., magneto- aco ustic waves within 
the chromosphere (as discussed in Schroder & Cuntz| 
l2QQ5f ). In contrast, models with larger stellar masses, 
at roughly M > 1.0 M©, would quickly become C-rich 
as a consequence of the third dredge- up; under these 
conditions the long-period pulsation and dust condensa- 
tion would become efficient enough to trigger a radiative 
wind, with consequent strong enhancement of the mass 
loss and quick termination of the AGB phase. This pre- 
diction is in agreen ient with the results of d etailed models 
of dust formation (jFerrarotti fc Gaiill2006f ). 

Finally, we remark that we apply the above-described 
mass- loss formulas solely during the TP-AGB. We rec- 
ognize that it would be important to explore the effect 
they have also for the RGB and early- AGB phases, but 
this is difficult in practice. It would require the calcula- 
tion of more extended grids of evolutionary tracks, prop- 
erly evaluating the mass loss during the evolution on the 
early- AGB, and including additional sets of horizontal 
branch models of smaller masses than presently available. 
This is very time consuming. Moreover, the evolution- 
ary effects expected would be much smaller than those 
described above. 

3.5. Results with the new TP-AGB tracks 

Figures [5] and [6] show the final result of using the new 
TP-AGB tracks in the simulation of ANGST galaxies, 
for cases A and B, respectively. It is evident that the 
reduction of TP-AGB lifetimes, and consequently of the 
TP-AGB termination luminosity, leads to a much better 
description of the data in both cases. For some galaxies 
such as ESO294-010, M81F6D1 and M81K64, the model- 
data agreement can be qualified as excellent, with simu- 
lated AGB numbers within the 67 % confidence level of 
Poisson fluctuations in the data, and a good description 
of the luminosity function for AGB stars. 

In a few other cases, however, the agreement is still not 
completely satisfactory, and the models tend to overes- 
timate the AGB numbers by factors of up to 2, as for 
DD0113, DD044, IKN, and DD078. Notice that in ah 
these cases the observed AGB/RGB ratio is very small 
and close to its minimum value, which suggests that 
these galaxies are really dominated by very old popu- 
lations. Therefore, we consider that the mismatch in the 
AGB/RGB ratios from the models might be attributable 
to the errors in the SFHs and metallicities of these galax- 
ies, added of course to the residual mismatches in the 
stellar models of higher mass, and in the prescription for 
dust obscuration. 



Finally, for the galaxies ESO0540-032, M81K61, SCL- 
DEl, DD071, and M81F12D1, there is in general just 
a modest excess (smaller than ~ 50 %) in the predicted 
numbers of AGB stars. It is remarkable that this excess 
is generally more evident for the bright section of the 
AGB sequence, whereas close to the TRGB the numbers 
of predicted and observed AGB stars are in quite good 
agreement. We interpret the mismatch for bright AGB 
stars as being likely attributable to the younger popula- 
tion in these galaxies - with the causes being either in 
a slight overestimation of the lifetimes of the more mas- 
sive TP-AGB models or in a modest excess in the SFH 
derived from ANGST data for ages < 3 Gyr. 

From those galaxies in which there is a good model- 
data agreement, we can conclude the following: They 
indicate that the TP-AGB lifetimes above the TRGB 
are of the order of 1.2 to 1.8 Myr (cf. Fig. [3]) for 
the low-mass low-metallicity AGB stars (M < 1 M©, 
[Fe/H] < —1.2) which are typically being sampled. The 
need for such a reduction in TP-AGB lifetimes, with re- 
spect to those in MG07 ni o dels, has already been indi- 
cated by iGullieuszik et al] (|2008l ) in their study of the 
old metal-poo r Leo II dS ph, and to a lesse r exte nt also 
by IHeld et al.l (|2010D and iMelbourne et al.l (|2010[ ) while 
studying metal-poor galaxies with somewhat younger 
SFHs. On the other hand, it is quite reassuring that the 
final core masses of our best-fitting AGB models are be- 
tween 0.51 and 0.55 Mq (again for low-mass stars, with 
say M < 1 M0), which is in good agreement with the ob- 
served masses of white dw arfs in globular clus ters of the 
MW (of 0.53 ±0.01 Mo, cf . iKalirai et 1III2009I) . This lat- 
ter indication also agrees with the low initial-final mass 
relation derived from the number counts of hot white 
dwarfs - mostly fr om the thin di sk - in GALEX wide 
area surveys (Bianc hi et all 120101 ). Thus, we find that 
the dramatic reduction in the lifetimes and final masses 
of our TP-AGB models are well in line with the indica- 
tions from independent data. 

Whereas we have tested two different prescriptions for 
the mass loss rate at the dust-driven phase, both turn 
out to provide essentially the same results, whith case A 
reproducing our data for old metal-poor galaxies just 
slightly better than case B. However, the differences be- 
tween these two prescriptions are expected to increase 
for stars of masses > 2 M©, which reach significantly 
higher luminosities. Therefore, they need to be tested in 
galaxies containing young populations, which will be the 
subject of forthcoming papers. 

4. CLOSING REMARKS 

In this work, we have derived strong constraints on 
the lifetimes of low-mass metal-poor AGB stars, via 
the study of their numbers and luminosity functions in 
galaxy regions with a predominantly-old and metal-poor 
SFH. The AGB lifetimes at magnitudes brighter than the 
TRGB turn out to be between 1.2 and 1.8 Myr for stellar 
masses of M < 1 M©, at metallicities [Fe/H] < -1.2. In- 
deed, TP-AGB evolutionary tracks with these lifetimes 
are able to reproduce quite well the low AGB/RGB num- 
ber ratios in a subsample of these galaxies, as well as their 
luminosity functions. In some other cases, the agreement 
is still not perfect, but the new models perform signifi- 
cantly better than the previous MG07 ones. Although we 
have tested just two possible alternative prescriptions for 
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Fig. 4. — Evolution during the TP- AGB phase of a model with Mjnit = 0.832 Mq (corresponding to a mass at the first thermal pulse 
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iSchroder Cuntzi (|2005|) formalism. Refer to the text for more explanation. 



mass loss, they both produce a quite similar reduction in 
the AGB lifetimes, as well as similar final masses for the 
resulting white dwarfs. 

Whereas the constraints we derive may apply to quite a 
limited range of stellar masses and metallicities, they also 
represent an important step towards a robust calibration 
of the lifetimes of TP- AGB stars at all masses and metal- 
licities. Indeed, with this single point in mass-metallicity 
space we can discard some of the already-proposed for- 
mulas for the mass loss on the AGB, for instance the 



"Bowen & Willson (Tm) and ' WillsonI (|2QQQ[ ) one as im- 
plemented in MG07. Moreover, we can use this result 
as a reference point to start a more thorough calibra- 
tion of the AGB lifetimes at regimes of increasing masses 
and metallicities, making use of additional galaxies with 
the presence of younger SFH and more metal-rich giant 
branches. Such an extension of the AGB calibration, 
however, cannot proceed via the analysis of optical data 
only, as performed in this paper. It requires at least the 
use of near-infrared photometry in order to sample the 
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Fig. 5. — The same as Fig. El but now using our best TP-AGB models, for case A mass-loss. 



most massive and evolved AGB stars, as well as those 
of higher metallicities. Indeed, we are presently collect- 
ing near-infrared data with the new HST/WFC3 camera 
(IR channel), for some dozens of ANGST galaxies with a 
well-measured SFH. Forthcoming papers will discuss the 
problem of improving AGB models in the light of these 
new data. 

The present TP-AGB tracks surely represent an im- 
provement in our sets of isochrones and tools to simu- 
late stellar populations. For this reason they are being 
included as an alternative to the MG07 models, in the 



CMEEE and TRILEGAlO web interfaces, which provide 
isochrones and simulated photometry of resolved stellar 
populations in a wide variety of filter systems. 
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Fig. 6. — The same as Fig. [21 but now using the TP-AGB models for case B mass-loss. 
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